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Abstract 

The process e'^e" — > K'^ K~ J /ip is observed for the first time via initial state radiation. The 
cross section of e'^e~ — > K~ J /ip for center-of-mass energies between threshold and 6.0 GeV is 
measured using 673 fb~^ of data collected with the Belle detector on and off the T(4S') resonance. 
No significant signal for y(4260) — > K'^ K~ J /ip is observed, and we determine B{Y{A2QQ) 
K+K~J/i;)T{Y{A260) e+e") < 1.2 eV/c^ at a 90% confidence level. We also find evidence for 



e ' e 



KgKgJ/ip in the same data sample. 



PACS numbers: 14.40.Gx, 13.25.Gv, 13.66.Bc 
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The study of charmonium states via initial state radiation (ISR) at the S-factories has 
proven to be very fruitful. In the process e~^e~ —>■ 77S'i?vr+7r~ J/?/', the BaBar Collaboration 
observed the 1^(4260) [l]. This structure was also observed by the CLEO ^ and Belle Collab- 



orations with the same technique; moreover, there is a broad structure near 4.05 GeV/c 



in the Belle data. In a subsequent search for the Y{A2Q0] 
process, BaBar found a structure at around 4.32 GeV/c^ 



in the e^e — > 7/5_Rvr+7r ip{2S) 
Si, while the Belle Collaboration 



observed two resonant structures at 4.36 GeV/c^ and 4.66 GeV/c^ jsl. Recently, CLEO 
collected 13.2 pb~^ of data at \/s = 4.26 GeV and investigated 16 decay modes with char- 

n 

monium or light hadrons [6||. The large e^e tt+tt J/ip cross section at this energy is 
confirmed. In addition, there is also evidence for K~^K~ J/ip (3.7cr) based on three events 
observed. Further investigation on the process e'^e' K^K'J/ip will shed light on the 
understanding of the y(4260) and the other vector charmonium states. 

In this Letter, we use a 673 fb~^ data sample collected near the T(4S') with the Belle 
detector j?! operating at the KEKB asymmetric-energy e^e~ (3.5 on 8 GeV) collider 8] 
to investigate the K^K~ J /if) final state produced via ISR. About 90% of the data were 
collected at the T(4S') resonance (^i = 10.58 GeV), and the rest were taken at a center-of- 
mass (CM) energy that is 60 MeV below the T(45') peak. 

We use PHOKHARA jol that was validated in previous analysis 3] to generate signal 
events. In the generator, one or two photons are allowed to be emitted before forming the 
resonance X, then X decays into K~ J /ip with J/ip decays into e^e~ or fi~^fi~. When 
generating the MC sample, the mass of the X is fixed to a certain value while the width is 
set to zero. In X ^ K~^K^ J/tp, a pure 5- wave between the K^K~ system and the J/ip, 
as well as between the and K~ is assumed. The invariant mass of the K^K~ system 
is generated according to phase space. To estimate the model uncertainty, we also generate 
events with K^K~ invariant mass distributed like m^+^- in il){2S) ti^tx^ J/ip decays 10|, 
i.e., ^^'^'^ oc Phase Space x {m\+^- — Amj^y. 

We select candidate events with criteria similar to those used for the analysis of e'^e~ — 

n 

^J/ip [3[. We require the number of charged tracks to be four with a zero net charge. 
For these tracks, the impact parameters perpendicular to and along the beam direction with 
respect to the interaction point are required to be less than 0.5 and 4 cm, respectively, 
and the transverse momentum is restricted to be higher than 0.1 GeV/c. For each charged 
track, information from different detector subsystems is combined to form a likelihood for 
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each particle species (z), Li Tracks with7^ir = ^^f^^ > 0.6, are identified as kaons wit^ 
an efficiency of about 92% for the tracks of interest; about 4% are misidentified tt tracks 11| 
Similar likelihood ratios are formed for electron and muon identification. For electrons from 
J/tp e~^e~, one track should have TZe > 0.95 while the other track has TZ^ > 0.05, 
this results in a very pure J/tp e^e~ sample with an efficiency of 90%; for muons from 
J/ip ^ fi'^fi~, at least one track is required to have T^^j > 0.95; in cases where one of 
the tracks has no muon identification (ID) information, the polar angles of the two muon 
tracks in the K~ fi~^ fi~ center-of-mass system are required to satisfy —0.7 < cos6'^ < 0.7 
based on a comparison between data and MC simulation. The efficiency for J/tp ^ fJ'^fJ'~ 
is 87%. Events with 7-conversions are removed by requiring TZ^ < 0.75 for the K'^K~ 
tracks. For the J/ip ^ e^e~ mode, 7-conversion events are further suppressed by requiring 
a K^K~ invariant mass greater than 1.05 GeV/c^; this also removes events with a </> signal 
in the final state. For the J/ip ^ fi~^fi~ mode, we require a K^K~ invariant mass outside 
a ±10 MeV/c^ interval around the (f) mass to remove events with a (f) signal in the final 
state, possibly produced via e+e~ 77*7* — > ■y(f)£^£~. The detection of the ISR photon 
is not required, instead, we identify ISR events by the requirement \M^^^\ < 1.0 (GeV/c^)^, 
where M^^^ is the square of the mass that is recoiling from the four charged tracks. The 
M^^^ requirement is tighter than that in our previous analyses 0, Q] as the M^^^ resolution 
has improved, due to the lower momenta of the particles in K^K~ J/ip as compared to the 
Tx^-K" J/ip and TX'^'K~%l}{2S) final states. 

Clear J/ip signals are observed in both decay modes. We define a, J/ip signal region as 
3.06 GeV/c^ < nii+i- < 3.14 GeV/c^ (the bremsstrahlung photons in the e^e~ final state 
are included, and the mass resolution is about 17 MeV/c^), and J/ip mass sidebands as 
2.91 GeV/c2 < m£+£- < 3.03 GeV/c^ or 3.17 GeV/c^ < m^+t- < 3.29 GeV/c^; the latter are 
three times as wide as the signal region. 

Figure [1] shows the K'^K" J/ip invariant mass |l^ distribution after the above selection, 
together with the background estimated from the J/ip mass sidebands. There is a broad 
enhancement around 4.4-5.5 GeV/c^. In addition^here are two events near y/s = 4.26 GeV, 
where CLEO observes three K^K~ J /ip events |6|. It is evident from the figure that the 
background estimated from the J/ip sidebands is low, which indicates that the background 
from non-J I'll) final states is small. The backgrounds not measured from the sidebands, such 
as XJ/ip, with X not being K^K~ , are found from MC simulation to be less than one event 
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and are neglected. 
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FIG. 1: The K^K^J/ip invariant mass distribution. The open histogram is the selected data 
while the shaded histogram shows the normalized sideband events. The solid (dashed) curve shows 
the best fit with two (one) Breit-Wigner functions and an incoherent background term, while the 
dash-dotted curve indicates a fit to the sideband background. 

The data points in Figs. EJ^a) and (b) show the M^^^ distribution (the requirement on it 
has been relaxed) and the polar angle distribution of the K^K~ J/ip system in the e~^e~ CM 
frame for the selected K~^K~ J/ip events with invariant mass between 4.4 and 5.2 GeV/c^. 
The data agree well with the MC simulation (shown as open histograms), indicating the 
existence of signals that are produced from I SR. 

We estimate the significance of the events between threshold and 6.0 GeV/c^ by calcu- 
lating the probability that the estimated number of background events in the normalized 
J/ip sidebands (12.3 ± 2.0) fluctuates to the number of observed events in the J/ip signal 
region (93) or more. It is found that the above probability is very small, corresponding to 
a statistical significance for the signal much larger than lOo". 

The e~^e~ — > K~^K~ J/tp cross section for each K~^K~ J/ip mass bin is computed with 



^obs _ ^bkg 



where nf'^, ''^^^^ and Ci are the number of events observed in data, the number of 
background events from a fit to the J/ip sideband events, the efficiency, and the effective 
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FIG. 2: M^^gj, distribution (a) and the polar angle distribution of the K~ J/ip system in the 
e^e~ CM frame (b) for the selected K~^K~J/ip events with invariant masses between 4.4 and 
5.2 GeV/c^. The points with error bars are data, the shaded histogram is the normalized J/^p 
sideband distribution, and the solid histograms are MC simulated events, normalized to the mea- 
sured cross section and integrated luminosity. The background from J/ip mass sidebands has been 
added to the MC simulation. 



luminosity 



13| in the i-th K^K J/ip mass bin, respectively [ij]. Due to the low statistics. 



we fit the background distribution with a second-order polynomial and take the fit number 
as the background in each K~ J / if) mass bin. The dilepton branching fraction, B{J/ip — >■ 
= 11.87%, is taken from Ref. [l^. The resulting cross sections are shown in Fig. [3l 



where t 

events 



le error bars indicate the combined statistical errors of the signal and the background 



16| . The cross section we measure in the 4.25-4.30 GeV bin is consistent with the 



direct measurement at ^/s = 4.26 GeV by the CLEO experiment Gj. 

There are a few sources of systematic errors for the cross section measurement. The 
particle ID uncertainty, measured using the same method as in Ref. [sl with pure track sam- 
ples, is 4%; the uncertainty in the tracking efficiency for tracks with angles and momenta 
characteristic of signal events is about 1%/track, and is additive; efficiency uncertainties as- 
sociated with the J/ip mass and M^^^ requirements are also determined from a study of the 
very pure e~^e^ "^(SS) ir^iT^J/ip event sample. In this study we find that the detection 
efficiency is lower than that inferred from the MC simulation by (2.5 ± 0.4)% relatively. A 
correction factor is applied to the final results and 0.4% is included in the systematic error. 
Belle measures the luminosity with a precision of 1.4%, and the uncertainty of the ISR 
photon radiator is 0.1% 13|. The main uncertainty in the PHOKHARA generator 9| is 
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FIG. 3: The measured e~^e~ K~^K~ J /ip cross section for CM energies between threshold and 
6.0 GeV (points with error bars). The errors are statistical only; a 10% systematic error that is 
common to all the data points is not included. Bins without points have a central value of zero. 
The dashed curve shows the energy-dependent selection efficiency with scale in right-hand side. 

due to the modelling of the K^K~ mass spectrum. Figure H] shows the K^K~ invariant 
mass versus K~^K~ J/tp invariant mass, as well as the projection on the K^K~ invariant 
mass for events in the J/ip signal region. The K^K~ invariant mass tends to be large 
and close to the phase space boundary, with an accumulation of events at 1.2 GeV/c^ and 
1.7 GeV/c^. Simulations with modified K^K~ invariant mass distributions yield efficien- 
cies that are higher by 2-5% for different K~^K~ J/ip masses. This is not corrected for in 
the analysis, but is taken as the systematic error (conservatively estimated as 5%) for all 
K^K~J/ip mass values. The angular distributions of the final state particles are compared 
with the MC generation, no evidence was found for non-S-wave components. Estimating the 
backgrounds using different J/ip mass sidebands results in a change of background events 
at the 0.18/50 MeV/c^ level, corresponding to an average of about a 6% change in the cross 
section. According to the MC simulation, the trigger efficiency for the final state is 98%, 
with an uncertainty that is smaller than 1%. From Ref. 15|], the uncertainty on the world 
averages for B{J/ip = B{J/ip e+e") + B{J/ip — > determined by linearly 

adding the errors for the e~^e~ and fi~^fi~ modes, is 1%. Finally, the MC statistical error 
on the efficiency is 1.5%. These errors are summarized in Table [H Assuming that all the 
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sources are independent and adding them in quadrature, we obtain a total systematic error 
on the cross section of 10%. 
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FIG. 4: Scatter plot of nix+K- versus rrij^+j^-j/^ (a) and the projection on the K^K~ invariant 
mass (b) for events in the J/ip signal region. The shaded histogram is the distribution of the 
normalized sideband events. 



TABLE L Systematic errors of the cross section measurement. They are common for all data 
points. 



Source 


Relative error (%) 


Particle ID 


4 


Tracking efficiency 


4 


J/il) mass and M"^^^ selection 


0.4 


Integrated luminosity 


1.4 


mx+K- distribution 


5 


Background estimation 


6 


Trigger efficiency 


1 


Branching fractions 


1 


MC statistics 


1.5 


Sum in quadrature 
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An unbinned maximum hkelihood fit is apphed to the K~^K~ J/ip mass spectrum in 
Fig. [H Here the theoretical shape is multiplied by the efficiency and effective luminosity, 
which are functions of the K~^K~ J/ip invariant mass. The Breit-Wigner function for a spin 
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one resonance decaying into final state / with mass (M), total width (Ftot) and partial width 

to e+e" (Fe+e-) is 

s {s - M2)2 + M^rL PS{M) ' ^ ' 

where B{R f) is the branching fraction of the resonance to final state /, and PS{^/s) is 

the three-body decay phase space factor for X K~^K~ J/ip. The MC-determined mass 
resolution varies from 3 MeV/c^ at ntK+K-j/^ = 4.3 GeV/c^ to 6.8 MeV/c^ at 5.4 GeV/c^. 
This is small compared to the widths of the resonances in our study and is neglected. 

We fit the K~^K~J/ip invariant mass spectrum with one Breit-Wigner plus a back- 
ground term. The latter is a second-order polynomial that is fit to the scaled sideband 
data. The dashed curve in Fig. [1] shows the fit results. The resonance parameters are 
M = 4430^1 MeV/c2, Ttot = 2541^ MeV/c^, B{R K+R- J/tP) -Te+e- = 1.9 ±0.3 eV/c^, 
where the errors are statistical only. Although the peak mass is close to the ip{AA15), Ftot is 



larger than its world average of 62 ± 20 MeV/c^ 15|]. To determine the goodness of fit, we 
bin the data so that the minimum expected number of events in a bin is at least seven and 
determine a x^/nd/ = 10.7/6, corresponding to a confidence level (CL.) of 10%. Adding 
a coherent y(4260) amplitude in the fit with mass and width fixed at the Belle measure- 
ment Q] yields an upper limit on i3(F(4260) fsT+Js:- J/^)r(y(4260) ^ e+e") < 1.2 eV/c^ 
at 90% CL. A fit using the above functions together with a coherent ^/^(4415) component 



with mass and width fixed at its world average [15| values improves the fit (solid curve in 
Fig. [1]) to /ndf = 4.2/4, C.L.=38%; the significance of the '?/^(4415) signal is found to be 
around 1.7a with a branching fraction for '?/'(4415) K~^K~ J/ip at the few per mille level; 
in this fit, the mass and width of the second Breit-Wigner become 4875 ± 132 MeV/c^ and 
630 ± 126 MeV/c2, respectively. 

We also search for e~^e~ — > KgKgJ/ip with the same data sample. All the selection 
criteria are the same as for K~ J /ip except that the selection of K^K~ is replaced by the 
selection of two K^s decaying into tt^tt" 17|. After selection, the invariant mass distribution 
of the K^KgJ/tp candidates is shown in Fig. [51 Three events (one J/tp e~^e~ and two 
J/ip — >■ /U."'"/i~) are observed between 4.4 and 5.2 GeV/c^ where a large K~^K~ J/ip signal is 
observed. In the higher mass region, the number of events in the J/ ip signal region is about 
the same as expected from the normalized sideband events. 

MC simulation yields an average selection efficiency of e = (0.50±0.04)% for rrij^oj^oj/jp G 
[4.4, 5.2] GeV/c^. Assuming that there is no background, we obtain the average cross section 
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FIG. 5: KgKgJ/ip invariant mass distribution of the selected e~^e 



KgKgJ/tJj candidates, blank 



histogram is for events in the J/ip signal region, while the shaded histogram is the normalized 
sideband backgrounds. 

for e~^e~ — > KgKgJ/^ of o^kOrOj/iP = I-S^li pb, where the error is statistical only. The 
average cross section for e~^e~ K^K~ J/ip in the same mass range is ^k+k-j/^ = 3.1 ± 
0.6 pb, where the error includes combined statistical and systematic uncertainties. From 
the above two cross sections, we obtain R = = O.Glni, in agreement with the 

expectation {R = 1/2) from isospin symmetry. 

In summary, the process e'^e~ —>■ K~ J / ip is observed and the cross section is measured 
for the CM energy between threshold and 6.0 GeV. There is one very broad structure; fits 
using either a single Breit-Wigner function, or a ?/;(4415) plus a second Breit-Wigner function 
yield resonance parameters that are very different from those of the excited ip states currently 
listed in Refs. jisl. fisl]. We observe two events near the y(4260) mass, with a cross section 
consistent with the CLEO measurement at s/s = 4.26 GeV within the large errors. We 
set an upper limit on i3(F(4260) K+R- J/ ij)T(Y (4260) ^ e+e~) < 1.2 eV/c^ at 90% 
C.L. 
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